X-ray diffraction of a disordered charge density wave 
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We study the X-ray diffraction spectrum produced by a collectively pinned charge density wave 
(CDW), for which one can expect a Bragg glass phase. The spectrum consists of two asymmetric 
divergent peaks. We compute the shape of the peaks, and discuss the experimental consequences. 



The statics and dynamics of disordered elastic objects 
govern the physics of a wide range of systems, either peri- 
odic, such as vortex flux lines Q and charge density waves 
(CDW) 0, or involving propagating interfaces, such as 
domain walls in magnetic || or ferroelectric Q systems, 
contact lines of liquid menisci on rough substrates || and 
propagation of cracks in solids Q . It was recently shown 
that periodic systems have unique properties, quite dif- 
ferent from the ones of the interfaces. If topological de- 
fects (i.e. dislocations etc.) in the crystal are excluded, 
displacements grow only logarithmically 0, El El j instead 
of the power-law growth as for interfaces. The positional 
order is only algebraically destroyed 0, leading to di- 
vergent Bragg peaks and a nearly perfect crystal state. 
Quite remarkably, it was shown that for weak disorder 
this solution is stable to the proliferation of topological 
defects, and thus that a thermodynamically stable phase 
having both glassy properties and quasi-long range po- 
sitional order exists |lfj. This phase, nicknamed Bragg 
glass, has prompted many further analytical and experi- 
mental studies (see e.g. [ill LL2j for reviews and further 
references). Although its existence can be tested indi- 
rectly by the consequences on the phase diagram of vor- 
tex flux lines, the most direct proof is to measure the 
predicted algebraic decay of the positional order. Such a 
measurement can be done by means of diffraction exper- 
iments, using either neutrons or X-rays on the crystal. 
Neutron diffraction experiments have recently provided 
unambiguous evidence ^jjj of the existence of the Bragg 
glass phase for vortex lattices. 

Another periodic system in which one can expect a 
Bragg glass to occur are charge density waves 2].' where 
the electronic density is spatially modulated. Disorder 
leads to the pinning of the CDW jl4| . In such sys- 
tems very high resolution X-rays experiments can be per- 
formed The resolution is in principle much higher 
than the one that can be achieved by neutrons for vor- 
tex lattices, consequently CDW systems should be prime 
candidates to check for the existence of a Bragg glass 
state. However, compared to the case of vortex lattices 
the interpretation of the spectrum is much more compli- 
cated for two main reasons: (i) the phase of the CDW is 
the object described by an elastic energy, whereas the X- 
rays probe the displacements of the atoms in the crystal 
lattice (essentially a cosine of the phase); (ii) since the 
impurities substitute some atoms of the crystal, the very 
presence of the impurities changes the X-ray spectrum. 



This generates non-trivial terms of interference between 
disorder and atomic displacements 0, 0. It is thus 
necessary to make a detailed theoretical analysis for the 
diffraction due to a pinned CDW. The study of the spec- 
trum has been carried out so far either for strong pinning 
or at high temperatures 0, 0, 0, 0] . 

In this paper we focus on the low temperature limit 
where a well formed CDW exists and on weak disorder, 
for which one expects to be in the Bragg glass regime. We 
show that the diffraction spectrum consists in two asym- 
metric peaks. In contrast to previous assumptions jlfij . 
we show that the asymmetry is present also in the weak 
pinning limit. The peaks are power-law divergent, with 
an anisotropy in shape. This form is consistent with the 
Bragg glass behavior 0] . The asymmetry is a subdomi- 
nant power-law too, with an exponent that we determine. 
We also briefly discuss the role of unscreened Coulomb 
interaction for the CDW on the diffraction spectrum. 

The general expression 0] for the total diffraction in- 
tensity in a crystal is given by 

= ^ E e-**'*-^ (fifje-^-^)) , (1) 

where u, is the atom displacement from the equilibrium 
position Rj — ja, with a indicating the lattice constant, 
fi the atomic scattering factor and (...) denotes the dou- 
ble average over the disorder and over the thermal fluc- 
tuations. As an example let us first consider the case of 
fixed atoms (tt, = 0). We obtain: 



I(q) = f 2 Y,S(q-K) 



AfNj, 



(2) 



K 



where A/ = // — / is the difference between the impurity 
/ and the host atom scattering factors, Nj = ni(l — nr), 
with nj is the impurity concentration and / is the aver- 
age scattering factor. The usual Bragg peaks, in corre- 
spondence to the reciprocal lattice vectors K, arise from 
the first term in (|2J), the second term is responsible for a 
background intensity, called Laue scattering, due to the 
disorder. 

In a second stage we take into account displacements 
of the atoms related to the presence of a CDW. To this 
purpose, we consider an electron density characterized by 
a sinusoidal modulation: 



p{x) = p cos(Qx + 4>{x)). 



(3) 
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4> is the phase of the charge density wave and Q = 2fcf, 
where kp is the Fermi wave vector. The associated 
Hamiltonian writes: 

H = J d d x°- (V0(:r)) 2 ± V I d d xK(x)p(x), (4) 

where d is the dimension of the space. The first term in 
the Hamiltonian (@J represents the elasticity. The elas- 
ticity is in fact anisotropic 20] along the Q-direction: 

ffei. = J dxd d - x y C j- {d x <t>f + f (dycji) 2 , (5) 

where x \\ Q, and c\ 3> C2. The compression along x cor- 
responds to an increase of electric charge density and thus 
pays the price of Coulomb repulsion, while distortions 
along the remaining d— 1 directions are much easier. We 
are led back to (0} by redefining the spatial variables x' = 
x/y/ci and y' = y/y/c2, with c = {cicir 1 )^ . The main ef- 
fect in the diffraction spectrum is thus to make the shape 
of the peaks anisotropic, but this will not change the 
overall divergence. The local, but anisotropic, elasticity 
© is valid beyond the distance at which the Coulomb in- 
teraction between various parts of the CDW is screened. 
If this length is very large, or if one want to examine 
short range regime one should keep the q-dependence in 
the elastic constants. This leads to a more complicated 
behavior that we will only briefly discuss here and will 
be examined in details elsewhere |2l|. The second term 
in the Hamiltonian |0J reflects the effect of the disorder 
on the electron density. The Gaussian random function 
describes the impurity distribution and is charac- 
terized by the correlator E(x)S(y) = Ni5(x — y), Vq is a 
positive constant which measures the impurity potential 
and finally the sign + (— ) is related to the repulsive (at- 
tractive) interaction between the electrons and the local 
impurity. In the following, we restrict our analysis to the 
repulsive case, po is absorbed in Vq and we define the 
disorder strength D = V 2 N I . 

A density modulation is accompanied by a lattice dis- 
tortion u given at low temperature by 

u{x) = ^ V cos(Q:r + 4>{x)). (6) 

We are interested in the behavior of the scattering in- 
tensity I(q) near a Bragg peak (q ~ K). Since \5q\ = 
\q — K\ <C K, we can take the continuum limit i — > x and 
we obtain from QJ: 




where f r = jjf d d re- lSqr and /V = / + A/a rf / 2 S(§). In 
(0 we have applied the standard decomposition in cen- 
ter of mass R and relative r coordinates (x = R + ^ and 
y = R — 0. The integration over R has already been 



performed because u vary slowly at the scale of the lat- 
tice spacing. Assuming that in the elastic approximation 
displacements remain small (ui « JJj), one can expand 
Q) as powers of Ku . Developing up to the second order 
we get 17] : 

I{q) = 7 d + / a + 7 trl pi.,with (8) 
h=fq 2 J{<\)<~)), 

4 = -zgA/a^ 2 /^ l^{- r -)u{ r -)-^H- r ^ , 
7 tripl . = -zgA/V^£(-^)£(^)M^) - • 

While the contribution represents the intensity due to 
the atomic displacements alone, the contributions i a and 
7tri P i. are generated by the coupling between the disor- 
der and the displacement. The presence of a CDW is 
signaled by the formation, around each Bragg peak, of 
two satellites at reciprocal vectors K ± Q. In absence 
of disorder (D = and <f> ~ const.) the displacement 
term has the form 7j = f 2 q 2u o J2k ^(<7 + 7C ± Q) and the 
other terms are vanishing: in this case the two satellites 
have the same intensity and the broadening is absent . To 
interpret the experimental findings 0, 0, 0] , in par- 
ticular to explain the measured strong asymmetry [2{| 
between the peaks at K + Q and at K — Q, we need to 
account for the effect of impurities. 

In the literature the term 7 a was evaluated by means of 
models 0, 0, which describe the pinning by impos- 
ing a constant value 0o to the phase in © in proximity 
of each impurity, and 7 tr i p i. was conjectured to be neg- 
ligible [13, |22] • In that approach, the observed satellite 
asymmetry is seen as a clear sign of the strong disorder; 
in fact, 4>q is not constant and, for sufficiently large do- 
mains, one should have 7 a cx cos(</>o) ~ 0- To go be- 
yond this phenomenological approach and also deal with 
the weak disorder limit, in which one expects the B ragg 
glass, we use a Gaussian variational approach 0, |23| . 
We first perform the average over the disorder using the 
standard replica techniques. The replicated Hamiltonian 
corresponding to Q| is 

H cS . = J ^ r ^^(V0 Q ) 2 -^cos(0 a (r)-0 6 (r)), 

a a,b 

(9) 

where T is the temperature and the sum over the n 
replica has to be considered in the limit n — > 0. We stress 
that, moving from the Hamiltonian (0J to its replicated 
version we also need to change the correlation functions 
containing explicitly the disorder: we have, for exam- 

pie, (E(-§)«(5)) - - 5 ^E„<P-(-5)«i(§)> off .. After 
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FIG. 1: Intensities of the different contributions to satellite 
peaks. The more divergent term, 7d, is symmetric. I a and 
hvipi. are antisymmetric, with 7 a 3> Ttripi.- 



some manipulations and using we obtain: 

Id = 7\ Vo / [e- tQr + c.c] C d (r) , (10) 

J r 

I a = -jAfquoVN ia d D f [e- iQr - c.c] C a (r) , 



where C* d (r) = (e^S)-^-*)))^ and C a (r) = 
Tn b (e 4 ^°^-'~^ 6 ^~5^) eff are the positional correla- 
tion functions controlling the behavior of each contri- 
bution. We notice that the intensity of the peaks at 
q = Q + K and q = K — Q is symmetric, as in the case 
of a pure system, for the displacement term Id, but it is 
antisymmetric for J a . The sum of these two terms leads 
to an asymmetry of the peaks. Fig. ^ show the behavior 
of the different contributions. 

Following the method used in [lfj for flux lines in pres- 
ence of weak disorder, we can calculate the various terms 
in 110|l . We look for the best trial Gaussian Hamilto- 
nian Hq — J G a b(q)4 > a(q)4 l b(~q) in replica space, which 
approximates ©■ Defining 



B ab (r) = ((Mr) ~ M0)) 2 )o 



(11) 



2T 



G(q) - G ab (q) cos qr 



where G is the diagonal element of G a b, and using the 
Gaussian approximation, the positional correlation func- 
tions become C a (r) = -^J2a,b e and C d( r ) = 

B{r) 

e 2 , where B is the diagonal element of B a b- Two 
general classes of solutions exist for this problem: while 
the first class preserves the permutation symmetry of the 
replica (RS), the second class (RSB) breaks the replica 
symmetry. It has been shown [lfj that the stable solu- 
tion for d > 2 corresponds to the RSB class, while the 
RS solution remains valid at short distance. Ch is simi- 
lar to the correlation calculated for flux lines |lfj and will 



be discussed later. To evaluate the contribution of the 
interference between disorder and displacement we fac- 
torize the antisymmetric term C a (r) = x(r)Cd(r). We 
first consider the RS approximation: 



X(r) = ^ 



1 - e 



-Tf G c (q) cos qr 



(12) 



where G c — is the connected part of G a b. In d = 3 
we estimate 



C™(r) 



2ir 2 



(13) 



The triplet term can be evaluated in an analogous 
way, but it gives non-zero contributions only consider- 
ing higher order harmonic terms in the electron den- 
sity. Equation © becomes: p(x) = po cos(n(Qx + 
(f>(x))), with n = 1,2. As we have already found 
for 7 a , we get an antisymmetric term with a pref- 
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i 

2nT- 



oc 

En 
a. 



Af quoNid D and a correlation C tr i p i 



b.c 



-i(0 C (i)-20a(i) + 0b(-i))\ 



at low temperature we finally obtain 
^tri P l.l7J ~ -j—e 2. 



In d — 3 and 



(14) 



It is interesting to evaluate, at this stage, the relative 
weight of the two antisymmetric terms in a satellite 
peak. We introduce the Fukuyama-Lee length (or Larkin- 
Ovchinikov length) 013 R a = (cV^) 1 / 4 ^ (for d = 3 
R a = c 2 /D) such that <j> varies on scale given by the 
length R a . The ratio of the two intensity peaks is: 



I-tripl. 



7r/ V 



(15) 



For weak disorder R a ^> a it follows that I a ^> Ttripl. an d 
we thus need only to consider J a and id- 

Since for d = 3 the RS solution is unstable, to obtain 
the correct physics one has to look for the RSB method. 
Within this scheme the off diagonal elements of 

Gab(q) are parameterized by G(q,v) where < v < 1 
and the solution is characterized by a variational break- 
point v c . The form of the symmetric part is given in 



Cd(r) 



(16) 



where 



6| 



^— measures the strength of thermal fluc- 

7rca ° 

tuations, and rj ~ 1 is the Bragg glass exponent in d = 3. 
At low temperature one has I ~ R a . The algebraic be- 
havior of is controlled by small v (v < v c ). Values 
of v above the breaking point (v > v c ) give the small 
distance contribution. Finally one finds v c = %4>tJ- 

To fully characterize the spectrum it still remains to 
evaluate \( r ) m the RSB scenario : 



X(r) = ~ [1- 



dve 



-T J (G(q)-G(q,v)) cos qr } 



(17) 
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FIG. 2: Ratio between the RSB and RS solutions for x( r )- 
At large distance this ratio tends towards a constant value b, 
where b ~ 0.018. This means that the RSB solution affects 
x( r ) on ty by a multiplicative factor. 



Restricting to the case d = 3, we write: 



G(q)-G(q,v) = - 
c 



2 

q 2 + l- 2 ' PJ v/va ~'(q* + (l) 2 y 

(18) 

By integrating l|18|l over q and with some manipulations, 
(fP7|l becomes: 



dt- 



X(r) = — 



I I' dzexpf-Sir 3 C - e ~ rt ' 1 



o 



t 



(19) 



The low temperature behavior (I ~ i? Q ) of this term 
is sketched in Fig. [21 As for the Replica Symmetry case, 
we have C a (r) oc ~e a - . We can now compare the two 
terms: 

,R a , 



/ d (if + Q) = f KW 



)'' 



(20) 



I a (K + Q) = ^njAf^Nuf^Kuo /(^)^. 

V R a J r r r 

After executing the d-dimensional Fourier Transforms, 
we conclude that both terms are divergent: in particu- 
lar, 1^ oc ]_ v and I a oc q d\-i ■ This effect, shown in 
Fig. is a clear sign of a quasi-long range positional or- 
dered phase. We have found that the peak at K + Q is 
smaller than the K — Q one, as the potential between the 
impurity and CDW is repulsive (we would have the op- 
posite asymmetry in case of an attractive potential). We 
observe that for an ideal infinite resolution experiment, 
the symmetric term would be dominant, since C&{r) de- 
cays to zero less rapidly than C a (r). However, if the 
divergence in (|20|l is cut by the finite resolution of the 
experiment both terms should be taken into account be- 
cause Id is quadratic in the small parameter Kuq whereas 
I a is only linear. 

The powerlaw lineshape is obtained for a short range 
elasticity. If the Coulomb interaction is unscreened, as 



might be the case in fully gapped systems such as the 
blue bronzes, the dispersion of c\ should be kept in JSJl. 
In that case c\ (q) ~ o il a 2 , which leads to peaks diverging 
even faster than (j20J ;21j . 

On the experimental side few detailed diffraction spec- 
tra are available at the moment. One case is doped blue 
bronzes where the lineshape corresponding to the CDW 
has been obtained after substraction of a Friedel oscilla- 
tion contribution [J^. The observed asymmetry of the 
peaks would be compatible with both strong and weak 
pinning. However given the short correlation length ex- 
tracted from the data, this particular experiment is most 
likely still in the strong pinning regime. It would thus 
be highly desirable to have more detailed analysis of the 
lineshapes either in this compound, for different impurity 
concentrations, or in less disordered systems, where one 
can expect a Bragg glass behavior. 

We thank J. -P. Poujet and S. Ravy for stimulating dis- 
cussions. 
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